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1. Introduction 



Abstract. We study several fundamental operators in harmonic analysis related to Bessel 
operators, including maximal operators related to heat and Poisson semigroups, Littlewood- 
Paley-Stein square functions, multipliers of Laplace transform type and Riesz transforms. We 
show that these are (vector- valued) Calderon-Zygmund operators in the sense of the associated 
space of homogeneous type, and hence their mapping properties follow from the general theory. 

o 
o 

Q. 

. In their seminal article [H] Muckenhoupt and Stein investigated in a systematic way har- 

CN , monic analysis associated with ultraspherical expansions and their continuous counterparts, 

Hankel transforms. That paper is considered as a starting point of an important development 
' connecting harmonic analysis and discrete and continuous orthogonal expansions. Later many 

authors contributed to the subject by studying various questions in different settings, includ- 
■ ing in particular expansions into classical orthogonal polynomials; see [31 Section 1] for sample 

^ • references. In the recent years one can observe an increasing interest in harmonic analysis of 

orthogonal expansions, as confirmed by the attention of many mathematicians and numerous 
papers. 

In this article we study several fundamental harmonic analysis operators in the n-dimensional 
00 ' setting related to Hankel transforms. This framework is connected with the Bessel operator 

m ' 

A--A-^ — — 

(N 



where A G [0, oo)" is a multi-index. The operator A;^ will play in our considerations a similar 
role to that of the Euclidean Laplacian in the classical setting. It is formally self-adjoint and 
nonnegative in L'^(Wl,dfix), where R" = (0,oo)"' and 

X ; dfix{x) = ^ • . . . • x^^" dx, X G Rl. 

The spectral decomposition of A;^, or rather its suitable self-adjoint extension, is given via the 
Hankel transform, see Section [2] for details. Moreover, A\ admits the decomposition Ax = 
Yli^i-^ij where Di = d/dxi are the usual partial derivatives, and D* are their formal adjoints 
in L^(]R" , d/UA). Thus Dj, i = l,...,n, are naturally associated derivatives with the Bessel 
operator. 
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The main objects of our study are the followmg operators related to Ax (see Section [2] for 
strict definitions): 

• maximal operators 

VF,^: ||exp(-tAA)/||^^(^,), P,^ : f ^ \\exi> { - t^/A^) f\\^^^^^^, 

• Littlewood-Paley-Stein type square functions 

atXr'- exp(-tAA)/||^,(j(|„|/2+fe)._l^^), 



• multipliers of Laplace transform type 



poo 

T^^ : f^Ax exp{-tAx)f m dt, 
Jo 

/•CO 

Jo 



• Riesz transforms 

We treat all these operators in a unified way, by means of the Calderon-Zygmund theory. Our 
main result, Theorem 12. H says that these are either scalar- valued, or can be viewed as vector- 
valued, Calderon-Zygmund operators in the sense of the triple (W^jdfix, I " I)) where | • | stands 
for the ordinary distance. According to the terminology of Coifman and Weiss [11], this triple 
forms a space of homogeneous type (this means, in particular, that the measure nx possesses 
the doubling property). Consequences are then delivered by the general theory. In particular, 
we conclude mapping properties in weighted spaces. 

Typically the main difficulty related to the Calderon-Zygmund approach is to show suitable 
kernel estimates. Inspired by earlier ideas used in certain Laguerre settings [HI EH 121], we 
present a convenient and transparent technique based on an integral representation of the Bessel 
heat kernel that emerges from Schlafli's Poisson type formula for the modified Bessel function of 
the first kind, see Section [3l It is remarkable that a similar method has been developed recently 
by Nowak and Sjogren |16j in the more complex setting of classical Jacobi expansions. 

Our present results constitute a continuation and extension of many earlier investigations in 
the Bessel setting. This concerns the fundamental paper [14] as well as more recent articles, see 
for instance [U [21 [6l [T] [3 |9l |TOl |23] . In all the mentioned cases, the one dimensional situation was 
considered. The study of the multi-dimensional Bessel setting has been undertaken only recently 
by Betancor, Castro and Curbelo [3lll|, and by methods not involving the Calderon-Zygmund 
theory. The results of this paper are also related to those of [3l H]. 

The paper is organized as follows. Section [2] contains the setup, strict definitions of the 
investigated operators and statements of the main results. In Section [3] we gather various 
preparatory facts and results needed to furnish the proof of the main theorem, which is then 
done in Section [H 

Throughout the paper we use a standard notation with essentially all symbols referring to 
the space of homogeneous type (M" , df-ix, | • |)- Thus C^(M" ) denotes the space of smooth and 
compactly supported functions in M". By {f,g)dfi^ we mean f{x)g{x)dfix{x) whenever the 
integral makes sense. Further, LP{wdfix) stands for the weighted space, w being a nonnegative 
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weight on M" . Given 1 < p < oo, p' is its adjoint exponent, 1/p + 1/p' = 1. For 1 < p < oo, 
we denote by = Ap{W^,dfix) the Muckenhoupt class of Ap weights related to the measure 
d/XA- While writing estimates, we will use the notation X < Y to indicate that X < CY 
with a positive constant C independent of significant quantities. We shall write X ~ y when 
simultaneously X <Y and Y < X. 

2. Preliminaries and main results 
Let A G [0, oo)". For z G M", consider the functions 

n 
i=l 



where J,y denotes the Bessel function of the first kind and order i/, cf. |25|. It is well known that 

|2 _ ^2 I I ^2 



for each z G M" , the function ip^ is an eigenfunction of the n-dimensional Bessel operator A; 
and the corresponding eigenvalue is |zp = + . . . + z^, 

Aav^^ = |z| V^, z G M!^. 
The n-dimensional Hankel transform hx defined by 

hxf{z)= [ ip^,{x)f{x)dfix{x), ze^l, 



plays in the Bessel context a similar role as the Fourier transform in the Euclidean setting. It is 
well known that hx is an isometry in Lp'idux) and it coincides there with its inverse, = hx- 
Moreover, for sufficiently regular functions /, say / G C^(M"), we have 

hx{^xf){z) = \z\''hx{f){z), z&Wl. 

Note that in dimension one and for A = one recovers here the setting of the cosine transform 
on the positive half-line. 

We consider the nonnegative self-adjoint extension of Aa (still denoted by the same symbol) 
defined by 

(1) Aa/ = hx{\z\^hxf). f G Dom(AA), 
on the domain 

Dom(AA) = {/ G L\df,x) : \z\^hxf G L^dfix)}- 
Then the spectral decomposition of Aa is given via the Hankel transform. 

The semigroup {M/^/} = {exp(— tAA)} generated by — Aa is usually referred to as the heat 
semigroup associated with the Bessel operator, or simply the Bessel heat semigroup. It has the 
integral representation 

(2) W^f{x)= [ Whx,y)f{y)diix{y), f e L^{dfix), xeRl, t > 0, 



where the associated heat kernel is given by (see [251 p. 395]) 
W,\x,y)= [ e-*l^lV.'(^)^^(y)dMA(^) 



(2^exp ( - ^(kp + 1,1^)) n(-^^r^'+^/^/A.-i/2(^), x,yeRl, t> 0; 
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here denotes the modified Bessel function of the first kind and order u (cf. [25j ) . It may be 
deduced that the integral in ([2]) actually converges for more general functions from weighted 
spaces with Muckenhoupt weights, producing always smooth functions of (x,t) G M" x ]R_|_, and 
thus providing a good definition of on these spaces, see Lemma 13.51 below. Moreover, ([2]) 
defines a symmetric diffusion semigroup in U'{dfj.x), I < p < oo, in the sense of Stein [221 P- 65], 
which is Markovian; see for instance |18l Section 6]. The semigroup {P/} = {exp{—t^/A\)} 
generated by the square root of the Bessel operator is called the Poisson-Bessel semigroup. By 
the subordination principle, it is related to {VF/"} by 

P^f{x)= W^^^ x^n^ t>0. 

Jo V™ 

The maximal operators of these semigroups are defined by 

W,^f = sup \Wt^f\, P^f = sup \P^f\. 

i>0 t>0 

As is well known, mapping properties of and P^ are connected with the boundary behavior 
of the semigroups. Notice that P^f < W^f, by subordination. 

Littlewood-Paley-Stein square functions based on {W^/} and {P^} have the general form 

where 2<r<cxD,A;GN, mG N*^ is a multi- index, |m| = mi + . . . + is its length, and 
9^ = (9™^ . . . 5™". Square functions of this form are important tools in harmonic analysis. 

Given a bounded measurable function M on M", the associated Hankel multiplier is 
defined by 

Tltf = hxiMhxf), feL^idfix). 

Clearly, is well defined in L'^{dii\). We say that M. is of Laplace transform type when it 
can be represented as 

poo 



M{z) = Mw{z) = \zY / e-l^l '^{s) ds, z G 



JO 

or as 

/>oo 

M{z) = Mp{z) = \z\ I e-l^l"V(s)rfs, z£Rl, 



for some bounded function ip on An important special case here is the choice ip{s) = s~^"' / 
r(l — 17), 7 real, producing the multiplier Aiw{z) = \z\'^'^'^ that corresponds to the imaginary 
power A^''^ of the Bessel operator. 

Riesz transforms in the Bessel setting are formally given, according to a general concept, by 

where m G N" and \m\ is the order of the transform. To give this definition a strict meaning, 
we introduce the space 

= {/ G C°°{Rl) : hxf G C^{Rl)}. 
This space is a dense linear subspace of Lp'{d^\). We have, see ([T|), 

A^'"i/v = /iA(i^ri'"|/iA/), f^c\ 
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If / G C^, then A^'"''/^ e C~(M!^) and therefore is well defined on C^. Then R^ can 
be extended uniquely to a bounded linear operator on Lp'{d^\). All these properties will be 
justified in detail in Section [4.41 
We shall treat all the operators 

Ti^A pA \W X,P rpX r>x 

in a unified way, by means of the general C alder on- Zygmund theory. Clearly, the maximal 
operators and the g-functions are not linear. They are, however, associated with vector-valued 
linear operators taking values in some Banach spaces B, where B = L''{t^''+\"'^/'^>-^dt) in case of 
9rnkr ^^'^ ^ ~ L*" in casc of (7^^ ^. For and we shall, for technical reasons, 

choose B not as L°°{dt) but as the closed and separable subspace Co C L^{dt) consisting of 
all continuous functions f on M_(_ which have finite limits as t — )• O"*" and vanish as t — )• oo. In 
all the cases we shall say that the operator is associated with the corresponding Banach space 
B. Similarly, the linear operators and i?,^ will be said to be associated with the Banach 
space B = C. To obtain mapping properties of our operators, we shall prove that they are 
C alder on- Zygmund operators in the sense that we now explain. 

Let B be a Banach space and let K{x,y) be a kernel defined on x IR"\{(x,y) : x = y} 
and taking values in B. We say that K{x, y) is a standard kernel in the sense of the space of 
homogeneous type {W^,dfix, | • |) if it satisfies the growth estimate 

""^^^'^^""^-,A(i5(x,l-,|)) 

and the smoothness estimates 

(4) \\K{x,y) - K{x',y)\\M<Y- — ^ /p/ ] m' |x - y| > 2|x - 

\x-y\ Hx[B[x,\x - y\)) 

(5) \\K{x,y)-K{x,y')\\M<^y^ \ ^, \x - y\ > 2\y - y'\- 

\x - y\ l^x{B{x, \x - y\)) 

here B{x,R) denotes the ball in M" centered at x and of radius R. When K{x,y) is scalar- 
valued, i.e. B = C, the difference conditions @ and ([5]) are implied by the more convenient 
gradient condition 

1 



(6) \^x,yK{x,y)\ _ \ fryf \ \W 

\x - y\fJ-xiB{x,\x - y\)) 

Suppose that, for some 1 < r < oo, T is a linear operator assigning to each / G L'''{d^\) 
a measurable B-valued function Tf on R" . Then T is said to be a (vector- valued) Calderon- 
Zygmund operator in the sense of the space (M+, d^\, | • |) associated with B if 

(a) T is bounded from U'{dnx) to U^{d^\), 

(b) there exists a standard B-valued kernel K[x, y) such that 

Tf{x)= K{x,y)f{y)dfixiy), a.e. x ^ supp/, 

for every / G LJf (M" ), where L^(M") is the subspace of ) of bounded measurable 

functions with compact supports. 

Here integration of B-valued functions is understood in Bochner's sense. For the theory of 
Bochner integrals we refer to [27]. The Bochner-Lebesgue spaces L^{d^\), 1 < p < oo, are 
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defined to consist of all strongly measurable functions /: M" — )• B such that < oo, 
where 

II .11 I l/lRn ||/(x)||^dA^A(2;)) , l<p<00 
ll/llLjJ(dMA) =1 

[esssup^gjjn ||/(x)||b, p = oo 

According to [12J and [14], we define the atomic Hardy space of B- valued functions on 

M" to consist of these / G L\(dix\) which admit the atomic decomposition 

oo 

the series being convergent in L^{d[x\)^ where G C are scalars such that \ < oo, and 

aj are atoms in the sense we now describe. A strongly measurable B-valued function a is an 
atom in the context of (M" , d[i\ , | • | ) if 

(Al) there exist xq € and tq > such that suppa C i?(xO)''o)) 
(A2) ||a||i^(iRn) < l/^A(S(a;o,ro)), 
(^3) a{x)diix(x) = 0. 

The norm in H^''^^ is defined by 

oo 

inf^lojl, 

i=i 

where the infimum is taken over all atomic decompositions / = Xlj^i '^j'^j °f /• Note that the 
condition (A2) above may be replaced by 

(A2') ||a||L^(iRn) < {fix{B{xo,ro))y/''-^ for some fixed 1 < r < oo. 

In the one-dimensional case and for B = C the atomic Hardy space in the Bessel setting was 
investigated in O [THl . 

The space BMO^{dfix) is defined to consist of all locally d/iA-integrable B-valued functions / 
on such that 

l|/llBMOB(dMA) = sup I \\f{x) - fsh dlMx{x) < OO, 

where the supremum is taken over all balls B in {W^,dfix, | • |), and Jb is the mean value of / 
over B, Jb = y^j^ Jb fiv) ^l^xiy)- When B satisfies certain mild conditions, the dual of -ffj^'^t 
is BMOB*(d^A), with B* being the dual of B. 

It is well known that a large part of the classical theory of Calderon-Zygmund operators 
remains valid, with appropriate adjustments, when the underlying space is of homogeneous type 
and the associated kernels are vector- valued, see for instance I2U] . In particular, if T is a 
Calderon-Zygmund operator in the sense of {W^,dfix,\ ■ |) associated with a Banach space B, 
then (see [8] Theorem 1.1] and references given there) 

(Ml) T extends to a bounded operator from LP{wdfix) to Ij^{wdfix)j for every 1 < p < oo and 
every w ^ Ap, 

(M2) T extends to a bounded operator from L^[wd^x) to weak L\{wdijLx)j for every w £ A\, 
(M3) T extends to a bounded operator from -f^^ to L\{d^x), 
(M4) T extends to a bounded operator from to 'QM.O-s,{d^x)- 
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In (M1)-(M4) it is implicitly assumed that T is given initially on dense subspaces being inter- 
sections of the relevant spaces with U'{djji\). 
The main result of the paper reads as follows. 

Theorem 2.1. Let A € [0,oo)", m G N", A: > 0, A; + |m| > 0, 2 < r < oo, = Mw or 
Ai = Aip. Then each of the operators 

TxrX pX X,W X,P rpX r>X 

can he interpreted as a C alder on- Zygmund operator in the sense of the space of homogeneous 
type (M!{:,o?//A,| • |) associated with a Banach space B, where B is Co, Co, U-{t'^^+\'^\/^>-^dt), 
U{t^^+\'^\>-^dt), C, C, respectively. 

Consequently, each of the operators listed above satisfies (M1)-(M4), with appropriately cho- 
sen B in each case. When B 7^ C, these mapping properties can be translated to the following. 

Corollary 2.2. Let A G [0,oo)", m G N", A; > 0, A; + |m| > and 2 < r < 00. Then 
the operators W^,P^,g^^^,g^^^^, viewed as scalar-valued operators, satisfy (M1)-(M4) with 
B = C. Moreover, the resulting extensions are given by the formulas defining initially these 
operators in L^{dfi\) (with r = 2 in case of the maximal operators), where the relevant integrals 
converge. 

Proof. The first statement is straightforward. The second one is justified by standard arguments, 
with the aid of Remark 13.61 below; see for instance the proofs of |17t Theorem 2.1] and \14\ 
Corollary 2.5]. □ 

The proof of Theorem 12.11 will be given in Section 21 Treatments of each of the operators are 
naturally divided into the following three steps. 
Step 1. The operator is bounded from L'^{dfi\) to L^{d^x) for some 1 < r < 00. 
Step 2. The operator is associated with an integral kernel in the sense of (b) above. 
Step 3. The kernel satisfies the standard estimates ([3]), ^ and 

In the next section we gather various facts and lemmas that will be needed in the proof of 
Theorem 12.11 

3. Preparatory facts and results 

The modified Bessel function Ii, has the following Poisson-type integral representation ob- 
tained by Schlafli, see [251 Chapter VI, Section 6-15]. For u > —1/2, 

(7) Iu{z) = e-'' dn,+,/2is), z > 0, 

where 0,^ is the measure on [—1, 1] given by the density 

(-[ _ s^^*?-! r/s 

' 0F2'?-i/2r(??) ' ' 



In the limit case ry = we put il.Q = ((5_i + <5i)/v27r, where 5_i and 5i are the point masses 
at —1 and 1, respectively, so that ([7]) holds also for u = —1/2. Then, for any A G [0,c«)"', the 
Bessel heat kernel can be written as 

(8) Wt\x,y) = ^^^^„/2+|A| / exp (- ^g(x,j/,g)) dn^is), x,y eU^, t > 0, 
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where |A| = Ai + . . . + A„, the function q is given by 

n 

q{x,y,s) = 1x1"^ + \y\^ + 2"^XjyjSj, x,yeRl, 

i=i 

and denotes the product measure (S)7=i on the cube [—1, 1]". 

The following result is a crucial point in our method of estimating kernels. It relates expres- 
sions involving certain integrals with respect to di}\ with the standard estimates for the space 
{M.l,dfxx,\-\)- 

Lemma 3.1 ([171 Proposition 5.9]). Assume that A € [0,oo)". Then 

^^,y,s))-/'-\'\-'/'dn,{s) < -. ^-^^ ^, x^y. 

\x-y\fix{B{x,\x-y\)) 
At this point it is perhaps interesting to observe that, see [17, Proposition 3.2], 

n 

Hx{B{x,R))^R''Y[{xj + Rf^\ xeRl, R>0. 

The result below will come into play when proving the smoothness estimates (jH) and ([5]) in 
cases when B 7^ C. It will enable us to reduce the difference conditions to certain gradient 
estimates, which are easier to verify. 

Lemma 3.2 ([24| Lemma 4.5]). If x,x',y E are such that \x — y\ > 2\x — x'\ and 6 = 

ax + {\ — a)x' for some a E [0, 1], then 

^q{x, y, s) < q{9, y, s) < 4q{x, y, s), s E [-1, 1]". 

The same holds after exchanging the roles of x and y. 

The following technical result will be used repeatedly while showing the relevant kernel esti- 
mates. To prove it, and also for future use, we introduce some additional notation. Given two 
multi-indices m,M E N", the relation M < m means that Mi < rrii for all i = 1, . . . ,n. We 
write [dxq)^ to denote the quantity {dx^q)"^^ • . . . • {dx^q)"^" ■ 

Lemma 3.3. Let ^4 > 0, m, r E N" and k ef^ be fixed. Then 

d^d^dl [t-^ exp ( - lg(x, y, s)) ) < t-^-^-i\m\+\r\)/2 ^ _ ^^(^^ ^ 

uniformly in x,y €z M", t > and s E [—1, 1]". 
Proof. Taking into account the inequality 

|a + 6s| < (a^ + 6^ + 2a6s)^/^ a,b>0, sE[-1,1], 
we see that it is sufficient to show that 

d!^dTdl(t-^e.p ( - ^g(x,2/,.))) = j;c,,M,i?,(.)t-^-'=-^-(H+M+|M|+l«l)/2 

X {dxq{x,y,s)y\dyq{x,y,s))^{q{x,y,s)y exp - ^q{x,y,. 
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where the (finite) summation runs over < j < k, AI < m, R < r, and Cj^M,R are polynomials. 
Proceeding inductively, we arrive at the formula 

(9) 5,Mt-^exp(-^g(x,y,s)H = aj,kt-^'^-^ {q{x,y, s)y exp ^^q{x,y, s) 

o<i<fc 

where aj^k £ On the other hand, given i,j = 1, . . . ,n, another inductive reasoning leads to 

(10) a,7a,>xp(-lg(x,y,.)) = Yl /3M.i?,(s.,^,)t-("^+'-^+^^+^^)/' 



0<Mi<mi 



X {dxM^^y^ s))'^^' {dy^q{x,y, s))^' exp - ^g(x,y,s)^, 



where PMi,Rj{si, sj) are polynomials of two variables (some of them being null polynomials). 
The conclusion follows by combining ([9]) with ()10p . □ 



The next result provides a useful decomposition of dj^ip^. This will be necessary to prove L^- 
boundedness of operators involving higher order 'horizontal' derivatives, that is the (/-functions 
and the Riesz transforms. 



Lemma 3.4. Let A € [0, cxd)". Given m G N", there exist numbers cj G M, j < m, such that 
(11) 9-v^^(z) = z"^Yl CjixzYip'^'-Hz), x,zG Ml, 

with the notation = z^'^ • . . . • z^" and {xzy = (xizi)''^ • . . . • {xnZny" ■ 

Proof. By the tensor product structure of ip'^ it is enough to prove the result in the one- 
dimensional case. Then we proceed by induction on m. Using the identity (cf. [25\ Chapter III, 
Section 3-2]) 

— (2;~''Ji,(z)) = -z^^Jj.+iiz), 

we get dx^xi^) = ~z{xz)iPx'^^{z)- Suppose that (fTTjl holds for certain m G N. Then, with the 
aid of the recurrence relation (cf. [25, Chapter III, Section 3-2]) 

Mz) = —{J^-i{z) + J^+i{z)) , u>0, 

it follows that 

m 



(IfL III \ 

= ^'"'-\E 2A + 2j-i M'"''"-^'''("^ + (xzy^'^-^^^Hz)) -Yc,{xzy^'p^^^^Hz)y 

Thus we conclude that (jlip is satisfied with m replaced by m + 1. □ 

To prove the remaining two lemmas we will need the estimate 
(12) kx{z)\<l, x,zeRl. 

This is easily obtained from the following basic asymptotics for the Bessel function Jj^ (cf. [25\ 
Chapter III, Section 3-1 (8), Chapter VII, Section 7-21]): for > —1 we have 

J,(z)~z^ z^0+, J,(z) = o(^) 



Z —5- oo. 
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Lemma 3.5. Let A G [0,00)" and assume that f G U'{wdijL\), w G A^, 1 <p < 00. Then the 
heat integral of f , 

W^f{x)= f WtHx,y)fiy)dfi^iy), x G M^, t > 0, 
is a well-defined C°° function of {x, t) G M" x Moreover, given m G N" and A; G N, we have 
d^d^W^f{x)= f d^d'lW^{x,y)f{y)d^ix{y), x G M!^, t > 0. 

Furthermore, if f £ L'^{dfix), then we also have 

d^d^,W^f{x)= f dTdHe-'^'^"diz))hxfiz)d^x{z), xGR^, t > 0. 

Proof. Let C ^ C M" and F C F C M+ be bounded subsets of and M+, respectively. 
Observe first that by ([8]) and the fact that q{x,y,s) > \x — yp, s G [—1, 1]"", we have 

^tHx,y)<^^^jl^^ew{-^^\x-y\^), x,yeRl, t>0. 
This implies the estimate 

Wt\x,y)<e-^\y\\ yGRl, x e E, t e F, 

where c > is a constant depending, in particular, on E and F. Moreover, combining ([8]) with 
Lemma 13.31 we see that, given m G N" and G N, 

(13) \dTd!^W,\x,y)\<e"'^y\\ ye^l, x e E, t G F. 

Next we observe that the function y 1— t- e""^'^'^ belongs to all LP{wdfi\), w G A^, 1 < p < 00. 
This can be easily verified by splitting M" into dyadic 'rings', 

00 

M!^ = {x G M!^ : \x\ <l}u\J{xeMl: 2^ < \x\ < 2-''+^} 

j=0 

and using the estimate 

w{B{0, R)) = [ w{x) dfixix) < R>1, 

see for instance \15\ Section 4]. Similarly one can show that the function y 1— )• e"'^^^^'^ /w{y) is 
essentially bounded if w; G and to do that one uses the dyadic decomposition above and the 
estimate 

IIXiJ(0,fl)^~'l|oo<ii"+'l^l, R>1, 

following from the A^ condition, see Section 4]. 

Now let / G LP{wdfix), w G ^4^, 1 < p < 00. By the above observations and Holder's 
inequality we can write 

/ W^{x,y)\f{y)\d^,x{y)< ( e-'^\y\'\f{y)\di,x{y) 

II I I 2 1 1 

< ll/llLfMA'A)ll'2" llLp'(«,-p'/Pd;,;,)' xeE, teF, 

and the last expression is finite since w""^ G A^, (for p = 1 the L/ norm here must be replaced 
by lloo)- Thus is well defined for functions from the weighted spaces. The 
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fact that Wt^f{x) is a smooth function of {x, t) S M" x and the possibihty of passing with 
the derivatives dJ^df under the integral sign are proved inductively, by considering one partial 
derivative at a time and then using ()13p and the dominated convergence theorem. Here, and 
also elsewhere, the possibility of applying this theorem for differentiating under integral signs is 
justified with the aid of the Mean Value Theorem. 

Finally, we verify the last statement of the lemma. For / € L'^{dfix) we have 

W,^fix) = hx{e-'\'^"hxf){x)= [ e-'\'\\^,{x)hxf{z)df,xiz), x eMl- 
Further, by Lemma 13.41 and ()12p 

for some constant c > 0, and the function z i— )• e"'^'^' h\f{z) is integrable against dfix, as can be 
easily seen by means of Holder's inequality and the L^-boundedness of h\. In this position we 
can proceed inductively as before. □ 

Remark 3.6. From the proof of Lemma \3. 5\ we can conclude immediately the following con- 
vergence result. Assume that A G [0,00)*^ and fn ^ f in U'{wdpL\) for some w G A!^ and 
1 <p<oo. Let m E N" and k be given. Then, with t > fixed, 

d^d^Wt^fN d^d'^W^f, as N ^ 00, 

pointwise, and even uniformly on bounded subsets E C E C M" . 

Lemma 3.7. Let X G [0,00)'' and assume that f € C^{Wl). Then hxf G C°°(M") and, given 
m e N", 

dThxf{z)= [ 9,>^(x)/(x)d/iA(x). 

Proof. The arguments are based on Lemma[331 the estimate ()12p and the dominated convergence 
theorem. We leave elementary details to the reader. □ 



4. Proof of Theorem 12.11 

In this section we prove our main result. Theorem 12.11 In the following subsections we treat 
separately the cases of the maximal operators, the (^-functions, the Laplace transform type 
multipliers and finally the Riesz transforms. 

In the sequel we will often omit the arguments and write shortly q instead of q{x,y,s). We 
will also make a frequent use, without further mentioning, of the fact that for A > and 
B > 0, supj>Q t'^e"^* = Ca,b < 00. We shall sometimes tacitly assume that passing with the 
differentiation in x, or y or t, under the integral against di^\ or against dt is legitimate; similarly 
for changing orders of integrals. This is indeed the case in all the relevant cases, which may be 
verified in a straightforward manner by means of the estimates obtained in the proof of Theorem 
12.11 and the dominated convergence theorem. On the other hand, in more subtle cases we will 
always comment in detail operations of this kind. 
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4.1. Maximal operators and P^. 

First, recall that {Wt^} is a symmetric diffusion semigroup in the sense of Stein's monograph 
[22] . Then the maximal theorem [22l p. 73] applies showing the boundedness of on U'{d^\), 
1 < p < DO. The same is true for P^, by subordination. Theorem 12.11 complements these results 
by admitting weights and providing some further mapping properties of the maximal operators, 
see Corollary 12.21 

We treat in detail only the heat semigroup maximal operator and then only indicate how to 
make the arguments go through in case of P^. 

Step 1. We want to view as a vector- valued operator W'*', assigning to any / G L'^{djjL\) 
the function 

and bounded from L'^{dfix) to L'^^{dii\). We first ensure that indeed takes its values in 
L'jj^{d^\). Here the arguments are analogous to those from the proof of |l7t Theorem 2.1]. To 
make them work one needs two ingredients. The first one is the fact that, given / G L^idjix) and 
X G the function Wif{x) is continuous in t G (0,oo), which we already know (see Lemma 
13. 5|) . To get the remaining ingredient, it is enough to check that for / G L?'{d^x) 

(14) lim W^fix) = fix), lim W^fix) = 0, a.e. x G M" . 

This, however, is rather straightforward. We have 

Wt^f = hx{e-'\'\\xf), fGL\d^x), t>0. 

Then ()14p for / G follows by the dominated convergence theorem. Since is dense in 
L'^{dfj,x) (see Section 23] below) and the maximal operator is bounded on L'^{d^\), standard 
arguments show that (fH|l holds for / G L?'{dii\). 

Thus is a linear mapping from L'^{d^\) to L'^^^{dfi\), and as such is bounded, by the 
corresponding property of the scalar- valued operator W^. 

Step 2. The fact that is associated with the vector- valued kernel {W^{x,y)}t>o is justified 
exactly in the same way as the corresponding fact in the proof of [171 Theorem 2.1]. 
Step 3. We prove the standard estimates for the kernel {W^'{x,y)}t>o- By ([8]) we have 



Wt\x,y)\\L^(^dt) = sup-— ^ y"exp ( - ^) d^xis) 



(2t)"/2 

£ / wtwS (')"""" w^^^^W 



and the growth estimate ([3]) with B = Co follows from Lemma l3.ll 

To show the smoothness estimates dH) and ([5]) it is enough, by symmetry reasons, to consider 
Then using the Mean Value Theorem we get 



W,Hx,y)-Wt\x',y)\ < \x - x'\ V^Wt\x,y 



t\1 



\x=e 



where is a convex combination of x and x' that depends also on t. Thus it suffices to show the 
estimates 

l|9x-VF'?'(a;, v)| all ror/j^A 5, i ; rrr, — ; i = l,...,n, 
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for all x,x',y satisfying \x — y\ > 2\x — x'\. Using ([8]) and applying Lemma [3^ we get 

dnx(s) 
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< 

fn/2+ 



Now the desired bound follows by Lemma 13.21 and Lemma l3. 11 This finishes Step 3. 

Treatment of goes along the same lines. One has to combine the arguments given above 
with the subordination principle. In Step 1 the relevant identity is P/"/ = h\{e^*^^^hxf), f G 
Lp'{dfi\). We leave further details to interested readers. 

4.2. Square functions g^^^ and g^^^y 

We analyse in detail only the square function g'!^^^ based on the heat semigroup. The Poisson 
semigroup based ^-function is treated in a similar way, by means of the subordination principle 
(see for instance [24, Section 4.3]), hence the details are omitted. 

?mTr ^ vector-valued operator 



Step 1. We interpret g^^^ as the vector-valued operator , assigning to any / G L^{dfix) 



the function 



t>0' 



We will show that is bounded from U'{d^\) to L^{dfj,x), 2 < r < oo, where B = 

j^r ^^{k+\m\/2)r-i Actually, this task amounts to showing that the gf-function g^^ ^ is bounded 
on U'{dfj.x), 2 < r < oo. 

We first deal with the case r = 2. Let / G Lp'{di-Lx)- By Lemma 13.51 and Lemma 13.41 we see 
that 



j<m 



dTd^W,'f{x) = Y.cj {xzyd+\z)z"^{-l)>'\z\''e''^^^ hxf{z)df,x{z). 



Then 



X,W / f \\\^ 



xzyip>^+^iz)z"'{-ir\z\-''e-'\'\^hxfiz)d^,x{z) 



kli2k-t\: 



t 



2fc+|m|-l 



dt dfix{x) 



< 



E 

j<m 

E 

j<m 



{xzyv^,+^{z)z^{-l)'^\z\''^e-'\^\\xf{z)dfixiz) 
"'-^z\'^''e-'^'^\xf)(.x)fdfix+j(.x)t^''+\"'\-^dt 



dfix{x)t^''+\"'\-^dt 



Since hx+j is an isometry on L'^{dfix+j), this implies 



9m,k,2U )\\L2(dfj.x) 



e-^t\-\'t^'^+\-^h^dt\hxf{x)\^d^ixix). 



The integral in t here is equal T{2k + |m|)/(2|xp)^'^+l™l and we conclude that 



5^l2(/)|l22(rf^^) < / jl-^j^^\hxf{x)\-'dfix{x) < WhxfWh^af.,) = \\f\\h(df.,y 



2,2m|^|4fc 



Notice that when there is no horizontal component, i.e. m = (0, . . . , 0), then ||5'^'^2(/)lli^(<iMA) 
is equal to ||/||L2(d^;^)7 up to a factor independent of /. 
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Considering the boundedness for 2 < r < oo, it is enough to show that it fohows from the 
L''-boundedness of g^^2 ^^'^ maximal operator W^. To begin with, observe that by ([8]) 
and Lemma [33] 



so t'^~^\"^\/'^\^^^^W^^{x,y)\ < W2t{x,y). Consequently, for suitable /, 

i^+H/2|a-a*^p^A_^(^)| <Vr2^J/|(^), ^^Mi, t>0. 

Using this observation and Holder's inequality we obtain, for / G U{d^\), 
II '^'^ ( f^W 

oo JJ. 



Jw\_ Jo ^ ^ t>0 ' 



X I / I supt'=+i"^i/2|a-afT^iV(x)|)' d/xA(x) ■ ' 



t>0 

^ II X,W ( f\\\'^ II T/T/-^! f I ir~^ 



Thus if g^^2 ^^"^ bounded on U{dfix), then so is 

Step 2. We verify the fact that is associated with the kernel {9™5j^VF/'(x, y)}t>o. Note 

that the integral connecting with the kernel must be understood in Bochner's sense, and 

the underlying Banach space is B = [t^^^\"^\/'^^^~^dt). By density arguments it suffices to 
show that 

(15) ({a™a,^W^,V},>o'^) = ( f {dTd^W,\x,y)}t>ofiy)dfixiy),H 

^ ' \ JR'l 

for every / G C7~(R!^) and H{x,t) = Hi{x)H2{t), where Hi G C^{Rl), H2 G Cc°°(M+) and 
supp/ n suppi^i = 0. Here (•,•) means the standard Banach space pairing between L^{dij,x) 
and its dual Ll4d^lx), with B* = L'-'(t(*=+l'"l/2)^-idt). 
We start by considering the left-hand side of (fT5]) . 



1 ./IB"! 



The second identity above follows from Lemma 13.51 The change of the order of integration in 
the first identity is justified by Fubini's theorem. Its application is indeed legitimate since by 
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Holder's inequality 



^ 1 1 5^l2 (/ ) 1 1 L2 (rf^^ ) 1 1 ^^1 1 1 L2 (d^.;, ) 



|^{fc+|m|/2){r-l)-l/2^, 



2|lL2(dt)' 



and this quantity is finite in view of the L -boundedness of ^ 2- 

Now we focus on the right-hand side of ()15p . Interchanging the order of integrals we get 



{dTd^W,\x, y)}t>of{y) dfixiy), H 







1 ./TO" 



which confirms that both sides of psp coincide. Here application of Fubini's theorem is possible 
since 



1 ./IB" 



< 



^^^^^W^\x,y)\\^d^ix{y)d^,x{x) 



I I 

J SUDD -ffl J SU 



supp //i J supp / 

(i^A(y) d^A(a:) < oo, 



supp Hi Jsupp / fJ-xiB{x, \x - y\)) 

where we made use of the growth estimate for the kernel {d!^d^W^{x,y)}t^Q proved in Step 3 
below and the fact that the supports of / and Hi are disjoint and bounded. 
Step 3. We verify the standard estimates for the kernel {d'^d^W^^{x,y)}t>o taking values in 
B = L''(*{'=+l'"l/2)'-idi), where 1 < r < oo. By ([8]) and LemmaEJwe have 

\dTd^W,\x,y)\ ' ^ 



< 



fn/2+\X\+k+\m\/2 

This combined with Minkowski's integral inequality gives 



exp 



8t 



d:^d^W,Hx,y) 



< 



1 



dn^is). 



dflxis) 



p/2+|A|+fc+|m|/2 

The change of variable 1 1— )• qt shows that the norm under the last integral is equal to q^"/^^!'^!, 
up to a multiplicative constant. Thus now the growth estimate Q follows from Lemma l3. 11 

To show the smoothness estimates, we focus ourselves only on ([5]); the reasoning proving (jl]) 
is essentially the same. By the Mean Value Theorem, 



d^d>^wHx,y) - dTd^W,\x,y')\ <\y- y'\ Vyd^d^W,\x,y) 



:5m, ofcTT/A 



y=e 



where is a convex combination of y and y' depending also on t. Thus is suffices to show the 
estimates 

1 



dy^d:^d^w,\x,y) 



y=e 



< 



\x - y\li{B{x, \x - y\)y 



i = l,...,n, 



for all x,y,y' satisfying \x — y\ > 2\y — y'\. Using ([8|) and Lemma [3^ we see that 



dy,d:^d^W,\x,y)\ 



< 



1 



^n/2+|A|+fc+|m|/2+l/2 



exp 



8t 



dQxis). 
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This together with Lemma 13.21 produces (under assumption |x — y| > 2\y — y'\ 



< 



^n/2+|A|+fc+| 



From here we proceed as in the proof of the growth estimate above. This leads to 



dy,d:^d^w,\x,y) 



y=e 



< 



1 



n/2+|A|+l/2 '^^a(s), 



and the desired bound follows by Lemma l3.1i Step 3 is complete. 



4.3. Multipliers of Laplace transform type T^. 

We will consider multipliers Ai of the form 

M{z) = Mw{z) = \z\ 



''ip{s) ds, z E 



where ip is a bounded function on M+. The arguments given below and the subordination 
principle allow to treat in a similar way the multipliers Ai = Aip obtained from the formula 
above by replacing |zp by \z\ (this change corresponds to replacing the heat kernel by the Poisson 
kernel in the expression defining the associated integral kernel). 

Step 1. The fact that the operator T^/ = h\{A4h\f) is bounded on L?{d^\) is clear, in view 
of the boundedness of M and Plancherel's theorem for the Hankel transform. 
Step 2. We will show that is associated with the kernel 



KU{x,y) 



i>{t)dtWt\x,y)dt, x,yeRl, x^y. 



By density arguments, it suffices to verify that 



(16) 



dux 



KMix,y)f{y)9{x) dfixiy) dfix{x) 



for all f,g& with disjoint supports. We focus on the left-hand side of (fT6]l . By the definition 

-.A 
M 



of Tj^ and Plancherel's theorem for the Hankel transform we get 



hx{Mhxf){z)g{z) dfix{z) 
M{z)hxf{z)hxg{z) d^xiz) 



(17) 



hxf{z)hxg{z)dfix{z)dt. 



To write the last equality we used Fubini's theorem. Its application is justified since 



m |z|^e-*l^l hxfiz)hxgiz) dfixiz) dt < 



\hxf{z)hxg{z)\ dfixiz) 



< 



\oo\\j\\L2(df,>,)\\9\\L^dfMx) < 

We next analyse the inner integral in ()17p . Plugging in the integrals defining hxf and hx'g and 
then using Fubini's theorem we arrive at 

\z\^e-'\'\'hxf{z)hxg{z)df,x{z) 
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|z|2e-*l^lV^(a:)v3^ (y) djjxiz) f{y)^d^ix{y) dfix{x) 



^ ./HPi ./IB" 



n Imn 



-dtWt ix,y)f{y)g{x) dnxiy)dnxix). 



The application of Pubini in the first identity above is legitimate since, taking into account (jl2|) 
and that t > 0, 



1 /IB" Imn 



z|2e-*l"l ^^,{x)ip^,{y)f{y)g{x) dfix{z) dfixiv) dfix{x) 



oo. 



The second identity above is a consequence of the equality 

-dt [ e-*l^lV,\x)(p^(2/)d^A(z)= / \z\^e-'\^\'ip^,ix)ip^,{y)df,xiz), t > 0. 



Here passing with dt under the integral can be easily justified by ()12p and the fact that for 
t > e > we have |zpe"*l^l^ < |zpe"'^l''l^ z £ M" . 
Summing up, we proved that 



^Mf^9)a^=-I ml I dtWt\x,y)f{y)gix)dfix{y)df^xix)dt. 



To see that this expression coincides with the right-hand side in (|16p it is now enough to inter- 
change the order of integrals. An application of Fubini's theorem is again possible because 



ImdtWtix, y)f{y)g{x) | dt dfixiv) d^ix{x) 



< ll/lloolblloo / / — rm~\ rrr dfixiy) dnx{x) <oo, 

Jsupp / Jsuppg f^x{B{x, \X - 2/1)) 

where we used the assumptions on / and g and the estimate 



obtained implicitly in Step 3 below. The verification of (jl6|) is finished. 

Step 3. We show the standard estimates for the (scalar- valued) kernel K'^{x,y). By ([8]) and 
Lemma 13.31 we have 



\dtWhx,y)\ < ^^^7^ /exp ( - A) d^xis). 
Using the boundedness of ip and changing the variable of integration 1 1— )• qt we get 



Now the growth estimate ([3|) with B = C follows by Lemma l3. 11 

To prove the gradient estimate ([6]) , by symmetry reasons we may consider only the derivatives 
dxi, i = 1, . . . ,n. Then Lemma 13.31 gives 

\d.AWt\x,y)\ < ^^^7^4^ / exp ( - A) dnxis) 



18 



J.J. BETANCOR, A.J. CASTRO, AND A. NOWAK 



and hence, proceeding as above, 

The conclusion follows by Lemma l3. 11 

4.4. Riesz transforms i?^. 

Recall that, see Section O we defined 

(18) Rif = 9-A-l"l/V = a"^/iA(|^rl"^l/iA/), / e C\ 

First of all, we ensure that this definition is correct and that is dense in L'^{d^\). Indeed, 
given / G C^^ hxf, and thus also |z|~I'"I/ia/, belong to C^{MJ^), by the definition of C'^. Then 
by Lemma 13.71 the function h\{\z\~^"^^hxf) is smooth, so the formula defining makes sense. 
For the density of C'^, observe that h\{C^ (MJ^)) C C'^; this fohows from the definition of C^, 
Lemma 13.71 and the fact that hx coincides with its inverse in L?'{dii\). Since C^(M") is dense 
in Lp'{dii\) and hx is an isometry there, the conclusion follows. 

Step 1. We verify that R^ is bounded from C L?'{d^\) to L?'{diJL\). In consequence, it 
extends uniquely to a bounded linear operator on L'^{dii\). 
Let / G C^. Using Lemma 13.71 and Lemma 13.41 we can write 



Rlf{x)= / ^^^^Az)\z\-\^\hxf{z)d^^y{z) 
(19) =^c,-x^-/' ^l+\z)\z\-\^\z"'+^hxf{z)d^,x{z) 



j<m 

Then, taking into account that the Hankel transform is an L^-isometry, we see that 

j<m 
jr<m 



Step 2. We show that R^ is associated with the kernel 



1 

^^n{x,y) = Y^^^j^J^ W^(x,y)tH/2-irfi^ x,yGM'^, x^y. 
By density arguments, it is enough to justify the identity 
(20) {Rif,g)^^^= [ [ Ri{x,y)f{y)'^df,x{y)df,x(.x) 



for all f,g £ C^(M" ) such that supp / n supp 5 = 0. 

First we focus on the left-hand side of (|2Up . Notice that here we cannot apply p8p directly 
to Rmf because / is not necessarily in C'^. To overcome this obstacle, we shall use a limiting 
argument to express {R^f, g)dfix by a sum of certain integrals. Let {/at} C be a sequence 
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approximating / in L?'{d^\). Clearly, in view of the L^-boundedness of i?^, {R^fN,g)d^^ 
{R^f, g)dfi^ as — )• oo. On the other hand, for each we can write (see (fT9|) ) 

{RifN,9),^=Y.c, I I {xzy^^,-^\z)z^\z\-\'^\hxfN{z)d^^^{z)W)d^ix{x). 
Then interchanging the order of integrals we get 

(21) . = E / z^+^\z\-\"'\hxfN{z)h^+A^~'9){z)dfsx{z). 

Here the application of Fubini's theorem is legitimate since 



\z-^+J\z\-\"^\h^f^{z)hx+j{x-^g){z)\dfixiz) 

< \\\z\-^'"^z'^hxfN\\L2^d^^J\^^+J(^''9)\\LHd^,+,) ^ WfNhHdt,,) \\g\\LHd^.,) < OO. 

This chain of estimates shows also that one can pass to the limit with in (j2ip . Thus we get 
{Rif,g),^^ = E ^^7" z"^+^\z\-\"^\h,fiz)h,+,{x-^g)iz)df,,iz). 

We next analyse the right-hand side of (j2Up . Taking into account the assumptions imposed 
on / and g and the estimate 

|W^(x,y)kH/2-irfi< 



HxiB{x, \x - y\)) 
proved in Step 3 below, we may apply Fubini's theorem to get 

Rt{x, y)f{y)g{x) d^ixiy) dnx{x) 



n .l-an 



1 



r(|m|/2) 



OO 



™l/2-M / e-'\^\\>^{z)^^Jz)df,x{z)f{y)g{x)df,x{y)dfixix)dt. 



We now focus on the three inner integrals entering the last expression. Observe that we may 
exchange dJ^ with the integral against d^\{z); this can be justified by means of Lemma [3^ (fT2]) 
and the dominated convergence theorem. Then an application of Lemma 13.41 leads to 



a- / e-*l^l ^^,{z)^'y{z)d,ix{z)f{y)g{x)d,jix{y)d,j.x{x) 

n ./lB>n ./lB>n 

'\^\\^ c,ixzy^^,+^iz)^^iz)df,x{z)fiy)'^df,xiy)d^ixix) 
V c,- / e~'\^\' z^+ihxf{z)hx+,{x-^g){z) df^xiz). 



j<m 

The last identity is a consequence of Fubini's theorem, and its application is possible because, 
in view of the assumptions imposed on / and g, ()12p and the fact that t > 0, 



e-tkl z"'+^x^^^+^z)^^iz)fiy)gix)\d^ix{z)d^^xiy)d^^x{x)<oo. 

1 ./IB" 
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Summing up, we proved that the right-hand side of ()20p is equal to 

1 /'OO p 

j<m ^l^l/"^) JO -f^l 

To finish the proof of (|20p it suffices now to justify the possibihty of exchanging the order of 
integrals in the last expression. This, however, follows readily by means of the estimate enabling 
the application of Fubini's theorem leading to ()2ip . 

Step 3. We prove the standard estimates for the kernel R^{x,y). By ([8|) and Lemma 13.31 we 
have 

KwHx,y)\ < ^„/,^|,Vi^i/2 /-P ( - I) ^^(«)- 

This implies 

and the right-hand side here was already estimated in the required way in Step 3 of Section 14.31 
above. 

To show the gradient bound, we observe that again by ([5]) and Lemma 13.31 
\V^,yd-W,\x,y)\ < ^,^,^1,1^1^1/,^,/, /exp ( - A) 

and consequently 

\V,,yRiix,y)\ < j ____exp ( - -) dtd^^s). 
Now the conclusion follows as in Step 3 of Section 14.31 
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